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A study of the formation of excitons as a problem of two Dirac particles confined in two-layer 
graphene sheets separated by a dielectric when gaps are opened and they interact via a Coulomb 
potential is presented. We propose to observe Bose-Einstein condensation and superfluidity of 
quasi-two-dimensional dipole excitons in double layer graphene in the presence of band gaps. The 
energy spectrum of the collective excitations, the sound spectrum, and the effective exciton mass 
are functions of the energy gaps, density and interlayer separation. The superfluid density Us and 
temperature of the Kosterlitz-Thouless phase transition Tc are decreasing functions of the energy 
gaps as well as the interlayer separation, and therefore, could be controlled by these parameters. 

I. INTRODUCTION 

The many-particle systems of the spatially-indirect dipole excitons in coupled quantum wells (CQW's) 
have been the subject of recent experimental investigations [ll-Q- These systems are of interest, in 
particular, in connection with the possibility of Bose-Einstein condensation (BEG) and superfluidity of 
dipole excitons or electron-hole pairs, which would manifest itself in the GQW's as persistent electrical 
currents in each well and also through coherent optical properties and Josephson phenomena 

Recent technological advances have allowed the production of graphene, which is a 2D honeycomb 
lattice of carbon atoms that form the basic planar structure in graphite [lo| . Graphene has been 
attracting a great deal of experimental and theoretical attention because of its unusual properties in 
its band structure [lT| - [l^ . It is a gapless semiconductor with massless electrons and holes which have 
been described as Dirac-fermions [l^. Due to the absence of a gap between the conduction and valence 
bands in graphene, the screening effects result in the absence of excitons in graphene. However, the gap 
in the electron spectrum in graphene can be opened by applying the magnetic field, which results in 
the formation of magnetoexcitons [l6j . The BEG and superfluidity of spatially- indirect magnetoexcitons 
with spatially separated electrons and holes in high magnetic field have been studied in graphene double 
layer [l^ and graphene superlattice [H, [l^ . The electron- hole pair condensation in the graphene-based 
bilayers have been studied in j20l - l23| . However, the effective mass of magnetoexcitons increases when the 
magnetic field increases and, therefore, the Kosterlitz-Thouless critical temperature of the superfluidity 
decreases with increasing magnetic field. 

In this paper we propose a new physical realization of an excitonic BEG and superfluidity in two parallel 
graphene layers, when one layer is filled by electrons, and the other one is filled by holes. We consider 
two parallel graphene layers separated by an insulating slab (e.g. Si02) and propose the formation of 
the excitons due to the gap opening in the electron and hole spectra in the two graphene layers. The 
advantage of the consideration of exciton formed by an electron and a hole from two different graphene 
layers, separated by an insulating slab, is that the dielectric slab creates the barrier for the electron- hole 
recombination which increases the life-time of the exciton compared to the exciton formed by an electron 
and a hole in a single graphene layer. 

There are different mechanisms of the band gap opening in graphene. Substrate-induced band gap 
opening in epitaxial grap hene is caused by the breaking of sublattice symmetry owing to the graphene 
substrate interaction [2J]. When graphene is epitaxially grown on SiG substrate, a gap of ~ 0.26 eV is 
produced [23 |. The electronic structure of graphene can be tuned by an organic molecule. The band gap 
can be opened in graphene due to the charge transfer between an organic molecule and graphene [l^ . It 
was demonstrated by angle- resolved photoemission spectroscopy that a tunable gap in quasi-free-standing 
monolayer graphene on Au can be induced by hydrogenation j26| . The size of the gap can be controlled 
by hydrogen loading and reaches '--^ 1.0 eV for a hydrogen coverage of 8% [1^. The band gap tuning in 
hydrogcnated graphene was also analyzed within the density functional theory [27| . 

The equilibrium system of local pairs of spatially separated electrons and holes can be created by 
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varying the chemical potential, using a bias voltage between two graphene layers or between two gates 
located near the corresponding graphene sheets. For simplicity, we also call these equilibrium local 
electron- hole pairs as indirect excitons. Excitons with spatially separated electrons and holes can be 
created also by laser pumping (far infrared in graphene) and by applying perpendicular electric field as 
for the CQW's We assume that the system is in a quasi-equilibrium state. Below we study the 

low-density regime for excitons, i.e. exciton radius a < n~^/^, where n is the 2D cxciton density. In this 
system the effective exciton mass can be controlled by the gap. The effective exciton mass can be small 
relative to the mass of free electron, and the Kosterlitz-Thouless transition temperature Tc controlled by 
the gap is expected to be the same order or relatively high compared to the coupled quantum wells case. 

Our paper is organized in the following way. In Sec. |ll] we present the Hamiltonian of the spatially 
separated electron and hole in two different parallel graphene sheets separated by a dielectric in the 
presence of the band gap. In Sec. IIIII we obtain the single-particle energy spectrum of dipole excitons in 
two-layers graphene and find the effective exciton mass. In Sec. IIVI we obtain the spectrum of collective 
excitation in the weakly-interacting gas of dipole excitons. The density of the superfluid component and 
the temperature of the phase transitions for the system of dipole excitons in two-layer graphene in the 
presence of a band gap are obtained in Sec. [V] Finally, the discussion of the results and conclusions 
follow in Sec. EB 



Let us consider a system of electrons and holes located in two different parallel graphene sheets. In this 
system electrons and holes move in two separate sheets with honeycomb lattice structure. We assume that 
excitons in this system are formed by the electrons located in the one graphene sheet and, correspondingly, 
the holes located in the other. Since the motion of the electron is restricted in one graphene sheet and 
the motion of the hole is restricted in the other graphene sheet, we replace the coordinate vectors of the 
electron and hole by their projections ri and r2 on plane of one of the graphene sheet. These new in- 
plane coordinates ri and r2 will be used everywhere below in our paper. Thus, we reduced the restricted 
3D two-body problem to the 2D two-body problem on the graphene plane. Each honeycomb lattice is 
characterized by the coordinates (r^, 1) on sublattice A and (r^, 2) on sublattice B with j = 1,2 referring 
to the two sheets. Then the two-particle wavefunction, describing two particles in different sheets, reads 
^'(ri, si; r2, S2), where ri and r2 represent the coordinates of the electron and hole, correspondingly, and 
si, S2 are sublattice indices. This wavefunction can also be understood as a four-component spinor, where 
the spinor components refer to the four possible values of the sublattice indices si,S2; 



In other words, the spinor components are from the same tight-binding wavefunction at different sites. 
Each graphene sheet has an energy gap. Obviously the gaps in these sheets are independent and in the 
general case we can introduce two non-equal gaps Si and ^2, corresponding to the first and the second 
graphene sheet, respectively. The gap parameters 61, 82 are the consequence of adatoms on the graphene 
sheets, which create a one-particle potential. 

The corresponding hopping matrix for two non- interacting particles, including the energy gaps 5i and 
82 on the first and second sheets, correspondingly, then reads 



II. EXCITON HAMILTONIAN 




(1) 




V 









5i + 52 



d2 

S1-S2J 



(2) 



In Eq. ^ di = hvpi—i]^ — ]^)^ '^2 = ^^^"("^aij ~ afj) ^^-^ the corresponding hcrmitian conjugates 
are d\ = hvF{—i-^^ + g^^); d\ = hvp{^i-^^ 4- where xi, yi and X2, 2/2 are the coordinates of 

vectors ri and r2, correspondingly, vp ~ \/3at/{2h) is the Fermi velocity of electrons in graphene, where 
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a = 2.566 A is a lattice constant and t « 2.71 eV is the overlap integral between the nearest carbon atoms 
psj . This Hamiltonian allows us to write the eigenvalue equation for two non-interacting particles as 

Ho 9q ^eo^o , (3) 

which leads to the following eigenenergies: 



eo{h,6i;k2,S2)=±^Jkf + k-l + 6f+Sl±2^{kf + df)m + S^)^±^kf+Sf±^k'l + 5l . (4) 

where fci and k2 are momentum of each particle, correspondingly. Eq. (|4]) gives the energy spectrum 
for two non- interacting particles in the presents of the non-equal gaps energies Si and 52- The energy 
dispersion is symmetrical with respect to the replacement of particles 1 and 2. When there are no gaps, 
(5i = and 62 = 0, as it follows from dU the energy dispersion is ±fci ± /c2- 

Let's now consider the electron and hole located in two graphene sheets with the interlayer separation 
D, and interacting via the Coulomb potential V{r) = —e'^/e\/r^ + where r is the projection of the 
distance between an electron and a hole on the plane parallel to the graphene sheets, e is the electron 
charge, and e is the dielectric constant of the dielectric between graphene sheets. Now the problem for 
the two interacting particles located in different graphene sheets with the broken sublattice symmetry in 
each sheet can be described by the Hamiltonian 



n = 



(-5i + 52 + V{r) d2 di \ 

4 -Si-62 + V{r) di 

dl Si+52 + V{r) d2 

V d\ 4 di~d2 + V{r)J 



(5) 



and the eigenvalue problem for Hamiltonian ([5]) is 

-H* = (6) 

where ^' are four-component cigenfunctions as given in Eq.([T|). 

The Hamiltonian ([S]) describes two interacting particles located in two graphene sheets and satisfies 
the following conditions: 

i) when the interaction between particles vanished V{r) = it describes two independent particles, 
each located in the separate graphene sheet, having two independent gaps energies related to the broken 
sublattice symmetry in each graphene sheet. 

ii) when the gaps in each graphene sheet vanish. Si = and S2 ^ the Hamiltonian describes two 
interacting particles in one graphene sheet [2^ (let us mention that for Si = S2 = and D = the 
Hamiltonian ([S]) is identical to the Hamiltonian (2) in Ref. [29j representing the two-particle problem in 
one graphene graphene sheet if the band gap is absent) if a two-body potential is e^/er or in two graphene 
sheets with the interlayer separation D, and interacting via the potential V{r) = —e^ /e^r'^ -\- . 

iii) when both gaps vanish i5i = and (52 = 0, as well as two-body potential V{r) = 0, the Hamilto- 
nian describes two non-interacting Dirac particles. It is important to mentioned that eigenenergies are 
symmetrical with respect of replacement particle 1 and 2. 

In Hamiltonian ([5|) the center-of-mass energy can not be separated from the relative motion even 
though the interaction V = V{r) depends only on the coordinate of the relative motion. This is caused 
by the chiral nature of Dirac electron in graphene. The similar conclusion was made for the two-particle 
problem in graphene in Ref. [l^ , where two particles in a single sheet were considered without gaps and 
D = Q. 



III. SINGLE EXCITON EIGENVALUE PROBLEM 

Since the electron-hole Coulomb interaction depends only on the relative coordinate, we introduce the 
new "center-of-mass" coordinates in the plane of a graphene sheet {x,y): 

R = ari + j3v2 , 
r = n - ra . (7) 
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Here the coefficients a and /? are to be determined later. Apparently we can use the analogy of the two- 
particle problem for gapped Dirac particles in two-layer graphene with the center-of-mass coordinates 
for the case of Schrodinger equation. The coefficients a and /? will be found below from the condition 
of the separation of the coordinates of the center-of-mass and relative motion in the Hamiltonian in the 
one-dimensional "scalar" equation determining the corresponding component of the wave function. 
We are looking for the solution of ([S]) in the form 



vI/,(R,r)=e^'C-%-(r) . (8) 



Let's introduce the following notations: 



/C-i- = /Ce'~ = ICx + ilCy , 
/C_ = /Ce ''^ = ICx — ilCy , 

e = tan-i(^^), (9) 

and rewrite the Hamiltonian ([S]) in a form of the 2x2 matrix as 

f 02 + Vir)ao-S,ao+S2CT3 d \ 

\ Ol O2 + V{r)ao - Siao + S^a^ J ' ^ ^' 

where Oi and O2 are given by 

Ol = hvp (a/C_ — idx — dy) tro = hvpalC-cro — hvpiidx + dy)ao, (11) 

,n f ^/C„ + idx + dy 

y piC+ + idx - Oy 

where x and y are the components of vector r, aj are the Pauli matrices, ctq is the 2x2 unit matrix, also 
dx = d/dx and dy = d/dy. Analysis of the operators (fTT|) and (fT2|) shows that the coordinates of the 
center-of-mass and relative motion can be separated. 

For <f>aa we can rewrite the eigenvalue problem as a one-dimensional equation (see Appendix [X]) : 

(Sr^ + Vir) - ) <t>aa = [e + h- 52] <t>aa - (13) 

The other components of ^ are given as: 

*6 = -(ecTQ - iD2 - 5iafi - (^acrs - V{r)afi)-^iD\^ a (14) 

for 



and 



i^r = ( ^-V^- )-(5^.-5.J-o (16) 



]=dx.ai+dy^a2 (17) 



with Pauli matrices Uj and 2x2 unit matrix uq. Moreover, we have 



1 "'"^ 



e + 5i + ^2 - V{r) + ^— ^(5^ + a^J 

e — Ol + ^2 



(i5a;2 - 9y.,)(j)aa ■ (18) 
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Assuming r <^ D and substituting the second-order series expansion for the interaction potential V(r) 
-Vq + into Eq. ([13]), where Vq = /{sD) and 7 = /{2eD^), we obtain 



+ 7*^ Uaa 



e + Si" 52 + Vq- 



(hvpIC) 
27" 



(19) 



The last equation can be rewritten in the form of the two-dimensional isotropic harmonic oscillator: 

(- J-i(e)V2 + 7r2) (baa = :Fo{e)4>aa , (20) 

where 

eihvF? 



2 {e^ ^{61+62)^) 



:Fo = e + S,-S2 + Vo-^-^. (21) 
The solution of Eq. pO|) is well known (see, for example, Ref. [Slj) and is given by 

^"^'^-2N,[^, (22) 




-^i(e) V -^iW 

where N = 2ni -I- ^2 -I- 1 with rii = 0, 1, 2,3,.. ., 712 = 0, ±1, ±2, ±3, . . . , ±ni are the quantum numbers 
of the 2D harmonic oscillator. 

After some straightforward but lengthy calculations (cf Appendix[B]) we obtain the following expression 
for the energy in quadratic order with respect to /C 

^ -.{hvplCf , (23) 

where /i = (5i + 82 and Ci = 2-yN'^{hvp)'^ . Thus, from (|23p we can conclude that the effective exciton 
mass M is given as a function of total energy gap Si + 62 and the parameter Ci oc as 

(24) 

The effective exciton mass M as a function of total energy gap + 82 and the interlayer separation D 
defined by Ea. (l24l) is plotted in Fig.[T] According to Fig.[l] the effective exciton mass M increases when 
the total energy gap 5i -\- 62 and the interlayer separation D increase. The three-dimensional Fig. ^) 
demonstrates dependence of the effective exciton mass on the total energy gap and interlayer separation. 
The dependence of the effective exciton mass M on the interlayer separation D is caused by the quasi- 
relativistic Dirac Hamiltonian of the gapped electrons and holes in graphene layers. Let us mention 
that for the excitons in CQWs the effective exciton mass does not depend on the interlayer separation, 
because the electrons and holes in CQWs are described by a Schrodingcr Hamiltonian, while excitons in 
two graphene layers are described by the Dirac- like Hamiltonian ([5]). 




IV. COLLECTIVE PROPERTIES OF DIPOLE EXCITONS IN A TWO-LAYER GRAPHENE 

After having found the mass and the energy for a single exciton in the separated double layer of 
graphene, we turn now to an ensemble of excitons in this structure. Due to interlayer separation D 
indirect excitons both in ground state (ni = 7^2 = 0) and in excited states have non-zero electrical dipole 
moments. We assume that indirect exciton interact as parallel dipolcs. This is valid when D is larger 
than the mean separation (r) between electron and hole along graphene layers D ^ (r) . 

The distinction between exciton and bosons manifests itself in exchange effects [1, [s^ - fs^ . These 
effects for exciton with spatially separated electron and hole in a dilute system na^ ^ 1 (n ^ D~^) 
are suppressed due to the negligible overlapping of wave functions of two exciton in the presence of the 
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FIG. 1: The effective exciton mass M in the units of free electron mass mo as a function of the total energy gap 
for the different graphene interlayer separations D (a), as the function of on interlayer separation D for different 
values of the total energy gap (b) and as the function of the total energy gap and graphene interlayer separation 
(c). 



potential barrier, associated with the dipole-dipole repulsion of an indirect exciton [8|. Two indirect 
exciton in a dilute system interact as U{R) = e'^D'^ /{eR^), where R is the distance between exciton 
dipoles along the graphene layers. A small tunneling parameter t due to this barrier is [l7j : 



t = exp 

where 




2M I -—] dR 



slog (s;i%V(27i-nM2e4D4)) 



is the characteristic value of the center-of-mass exciton momentum hJC defined as k = \/2AI p^^x, where 
^ex is the chemical potential of the system (see below). In Eq. (P5|) . rg = [2Me'^D'^ / k'^Y/'^ is the 
classical turning point for the dipole-dipole interaction, s = 4 is the spin degeneracy factor for the 
excitons and M is the effective exciton mass in the ground state given by Eq. ([M)) . Then the small 
tunneling parameter t has the form t ^ eyip[—2h~^{MY/^eDa~^^^]. Therefore, we neglect the overlap of 
the exciton wavefunctions in the limit of large layer separation D and consider the gas of excitons as a 
Bose gas. Consequently, at sufficiently low temperatures the dilute gas of excitons forms a Bose-Einstein 
condensate (ssl . [36j . Formally, the exciton gas can be treated by the conventional diagram technique 
for a boson system. In particular, the effective interaction of the dilute two-dimensional exciton gas (at 
na^ ^ 1) can be described by a summation of ladder diagrams (35| . From the latter we obtain an integral 
equation for vertex function T, depending on three momenta p, p', P and the frequency O, as 

r(p,p';P,n).t/(p-pO + ./J^|<£^^3)^i^j^ (y^+o). (25) 

where C/(p — p') is a dipole-dipole interaction in momentum representation. This equation is also repre- 
sented by diagrams in Fig. [5] The chemical potential of the system is given by 

^2 

Me. = ^ = nor(0, 0; 0, 0) = noFo . (26) 

Equation can be solved easily when the excitons occupy the ground state ni = 77.2 = 0. Then the 
energy spectrum of the exciton is given by e{P) = P'^ /{2M), where the mass M is given by Eq. (|24p . 

The specific feature of a two-dimensional Bose system is connected with the logarithmic divergence of 
two-dimensional scattering amplitude at zero energy 0, [sj, [s^] . A simple analytical solution of Eq. 
for the chemical potential can be obtained if nMe^D^ /{fV^e) ^ 1, which gives for the chemical potential 



2M sM log [sh^e^ / {2TrnAPe^D^)] 
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FIG. 2: The equation for the vertex F in momentum representation 

The solution of Eq. ([25]) at small momenta provides the sound spectrum of collective excitations 
e(P) = CgP with the sound velocity Cs — y^Ton/ {AsM) = ^ f^ex/M. The appearance of a sound 
spectrum is a consequence of the dipole-dipole repulsion. This sound spectrum of the collective excitations 
reflects the possibility for the existence of an excitonic superfluidity at low temperatures in a double layer 
graphene, provided that the sound spectrum satisfies to the Landau criterion of superfluidity [ssl. l36j. 



V. SUPERFLUIDITY OF DIPOLE EXCITONS IN DOUBLE LAYER GRAPHENE 

The dilute exciton gas which was discussed the previous section, consisting of electron-hole pairs on 
the graphene double layer, forms a collective state whose excitations are sound-like modes. This might 
be true at low temperatures, whereas at higher temperatures phase fluctuations can destroy this 2D 
collective state by creating vortex-like excitations (i.e. by unbinding vortex-antivortex pairs). The latter 
have short-range correlations which prevent a superfluid state. Therefore, superfluidity is only possible for 
temperatures below a critical temperature Tc- This critical temperature describes a Kosterlitz-Thouless 
transition Q and is defined as 

~ 2kBM ' ^^^^ 

where Us (T) is the superfluid density of the exciton system at the temperature T, and ks is Boltzmann 
constant. 

The function ns{T) in ()28|) can be found from the relation Us = n — n„, where n is the total density and 
Tin is the normal component density. We determine the normal component density following the usual 
procedure [35| . Suppose that the exciton system moves with a velocity u. At nonzero temperatures T 
dissipating quasiparticles will appear in this system. Since their density is small at low temperatures, 
one can assume that the gas of quasiparticles is an ideal Bose gas. To calculate the superfluid component 
density we find the total current of quasiparticles in a frame in which the superfluid component is at rest. 
Then we obtain the mean total current of 2D excitons in the coordinate system, moving with a velocity 
u: 

where / (e(P)) = (cxp [e(P)/(fcBr)] — 1)~^ is the Bose-Einstein distribution function. Expanding the 
expression inside the integral and leaving the first order by Pu/(fcBr). we have: 

^-^^ ^ 2mJ (2^^ Fe ~ 27rft2 Mcf " ' ^^"^ 

where C(f ) is the Rieniann zcta function (C(3) — 1.202). Then we define the normal component density 
Tiji as 1351 1 

(J) = n„u . (31) 

Comparing Eqs. pip and pop . we obtain the expression for the normal density n„, which implies for the 
superfluid density 

3C(3)fc|T3 

' ~ 2^S2 cW ■ ^^^^ 
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It should be noticed that the expression for the superfluid density Ug of the dilute exciton gas in the 
double layer graphene in the presence of the band gaps differs from the corresponding expression in 
semiconductor coupled quantum wells (compare with Refs. 0, [s^l by replacing the total exciton mass 
M = nie + rrih with the effective exciton mass M given by Eq. ((24)) ). 

Using Eq. ([32]) for the density rig of the superfluid component, we obtain an equation for the Kosterlitz- 
Thouless transition temperature Tc. Its solution is 



Tr = 




32 /sMksT^ 
27 V TTh?n 



1/3 



1- 1 



7^0 

c 

21/3 



- (33) 



where is the temperature at which the superfluid density vanishes in the mean-field approximation 
(i.e., n,(TO)=0), 



T = — — 



1 



6<(3) 



(34) 



The behavior of the Kosterlitz-Thouless transition temperature Tc as a function of the total energy 
gap, exciton concentration n and the interlayer separation D is presented in Fig. |31 using Eqs. (P^]) and 
P4p . As we can see in Fig. [31 Tj, increases when the exciton concentration n and increases and decreases 
when total energy gap and and interlayer separation increaseincrcascs. 




FIG. 3: Kosterlitz-Thouless transition temperature Tc as a function of the total energy gap Si +62, exciton 
concentration n and the separation between the two graphene layers: (a) demonstrates how the transition 
temperature depends on the total energy gap 5i + ^2 for the four different values of exciton mass concentra- 
tion : n = 3.0 ■ 10" cm"^ 4.0 ■ lO" cm"^; 5.0 • lO" cm"^; 6.0 • lO" cm"^; (b) shows the transition temper- 
ature dependence on the exciton concentration for ten different values of the total enertgy gap: {Si +62 = 
0.70eV; 0.75eV; O.SOeV; 0.85eV^); (c) exhibits the dependence of the transition temperature on the interlayer 
separation D for the different values of the total energy gap values: {Si + 82 = 0.60el^; O.TOeV^; O.SOeV; O.QOeV). 



VI. DISCUSSION 



We have considered an electron-hole pair with attractive Coulomb interaction, where the electron and 
the hole live in two different graphene layers separated by dielectric. The distance between these layers 
is tunable such that we can vary the strength of the Coulomb interaction. Moreover, we assume a band 
gap in the dispersion of the electron and the hole which is caused, for instance, by doping the graphene 
layers with non-carbon atoms. The electron-hole pair forms an exciton whose mass M depends on the 
sum of the two gaps and on the layer distance D, as given by Eq. ([M)) . This result is generalized to 
a dilute gas of such excitons, which experiences a repulsive dipole-dipole interaction. The latter does 
not pose a problem because the dipoles are fixed by the double layers and can only interact as parallel 
dipoles. This allows us to consider the dilute excitons as point-like bosons which can be treated in a 
conventional self-consistent approach for bosons, leading to an effective interaction which is defined by 
the integral equation ([25]) . A solution of the latter for point-like particles of mass M provides us a 
sound-like spectrum of the quasiparticles, which represents superfluidity. The advantage of observing the 
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exciton superfluidity and BEC in graphene in comparison with these in CQW's is based on the fact that 
the exciton superfluidity and BEC in graphene can be controUed by the gaps which depend on doping. 
Note that we considered the superfluidity in two cases: first, an equiUbrium system of electrons and 
holes created by the gates, and the second case is the electrons and holes created by the laser pumping 
such that the excitons are in the quasi-equilibrium thermodynamical state. A temperature Tc which is 
the critical temperature of a Kosterlitz-Thouless transition was obtained. There is a superfluid state for 
T < Tc and a normal state for T > T^- The value of this critical temperature is given by Eq. 



Using the value of the total energy gap ^ 0.26 eV from Ref. [2J| and a interlayer distance D ~ \Q nm 
we obtain for the critical temperature Tc ~ 0.1 K for exciton concentration n = 5 x 10^^ cm^^, while 
for a interlayer distance Z? = 3 nm the critical temperature « 1.3 K and for D = 1 nm the critical 
temperature becomes Tc ~ 7.5 K. 

The superfluid state at T < Tc can manifest itself in the existence of persistent ( "superconducting" ) 
electric currents with opposite directions in the graphene layers. The interlayer tunneling in an equilibrium 
spatially separated electron-hole system leads to interesting Josephson phenomena in the system: to 
a transverse Josephson current, inhomogeneous (many sine-Gordon soliton) longitudinal currents, [39| 
diamagnetism for the case of magnetic field B parallel to the junction (when B is less than a certain 
critical value Bd, depending on the tunneling coefficient), and a mixed state with Josephson vortices for 
B > Bel- In addition, taking tunneling into account leads to the order parameter symmetry breaking 
and to a change of the phase transition type. The interlayer resistance relating to the drag of electrons 
and holes can also be a sensitive indicator of the transition to the superfluid state of the electron hole 
system [13, El. The existence of a local superfluid density below Tc can be detected, for example, by 
measuring the characteristic temperature dependence of the exciton diffusion on intermediate distances 

The advantage of observing the exciton superfluidity and BEC in graphene in comparison with these 
in CQWs is based on the fact that the exciton superfluidity and BEC in graphene can be controlled 
by the gaps which depend on doping. Note that we considered the superfluidity in two cases: first, an 
equilibrium system of electrons and holes created by the gates, and the second case is the electrons and 
holes created by the laser pumping such that the excitons are in the quasi-equilibrium thermodynamical 
state. Another advantage is that graphene is much cleaner than typical semiconductors used for CQW's, 
where the roughness of QWs boundaries is crucial. Therefore, disorder is much less of a problem in double 
layer graphene. 

In conclusion, we propose a physical realization to observe Bose-Einstein condensation and superfluidity 
of quasi-two-dimensional dipole excitons in two-layer graphene in the presence of band gaps. The effective 
exciton mass is calculated as a function of the electron and the hole energy gaps in the graphene layers, 
density and interlayer separation. We demonstrate the increasing effective exciton mass with the increase 
of the gaps and interlayer separation. The dependence of the exciton mass on the electron-hole Coulomb 
attraction and interlayer distance comes from the Dirac-like spectrum of electrons and holes. Wc show 
that the superfluid density Ug and the Kosterlitz-Thouless temperature Tc increases with increasing 
excitonic density n and decreases with the rise of the gaps 6i and 62, as well as the interlayer separation 
D, and therefore, could be controlled by these parameters. As we mentioned before, the energy gap 
parameters Si and 62 are determined by the doping concentration. 

Appendix A: Eigenvalue Problem for two particles 

For the Hamiltonian PT7|) the eigenvalue problem "H^* = e^* results in the following equations: 
(©2 + Vir)ao - 6i(7o + 62(73) + d*,, = ecrn^a 
Ol-^a + {O2 + V{r)ao - Siao + (^afTs) = ecro^'fc ■ (Al) 

From Eq. (jAl[) we have: 

= {eao -02~ V{r)a^ + S^a^ - -Js^ts)"' 0\^a ■ (A2) 

Assuming the interaction potential and both relative and ccnter-of-mass kinetic energies are small com- 
pared to the gaps and 82 we use the following approximation: 

(eao -O2- V{r)ao + S^a^ - (^a^Ts)"' ^ -j^ j . (A3) 
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Using the fact that the operator 0\0i is purely herniitian, applying Eq. (|Aip and 



we obtain: 



(O2 + V{r)ao - Siao + (^ada) + fi^v^ 



2ia{K.xdx + JCydy)) 



ea-Q + Sictq - 620-3 



Now we rewrite Eq. (jASP in the following form: 



2„.2 



-5i + 62 + Vir) + h'v 



- V2 - 2ihvFa{ICa;d^ + JCydy) 



hvF il3IC^ + idx + dy) (t>ab = e(t>aa , 

hVF (/3/C+ + idx - dy) (t)aa + 



-Si -62 + V{r) + h^vl 



a'^K? - - 2ia{]C^d^ + JCydy) 



e- Si + S2 
We solve Eq. (jA7|) with respect to i/'ab^ 

a^JC^ - V2 - 2ia{JCA + JCydy) 



Oab = e^ab 



(A4) 
(A5) 

(A6) 

(A7) 
(A8) 



e + Si+S2- V{r) - h^v 



2„,2 



e - Si + S2 

Substituting ipab from Eq. (|A9|) into Eq. (|A6|) . we obtain: 



{I3JC+ + id., - dy) hvF^aa ■ (A9) 



-Si + S2 + V{r) + h^v 



2 2 a^lC^ - - 2m(/C:,9^ + JCydy) 



+ h'vj, {pJC^ + id, + dy) 
X {/3JC+ + id, - 9y) = e(^QQ 



e - (5i - (52 

e + oi + (52 - y(r) - h^vp ^ 



e — Si + S2 



(AlO) 



Assuming again that the interaction potential and both relative and center-of-mass kinetic energies are 
small compared to the gaps Si and S2 wc apply to Eq. (|A10p the following approximation: 



e + (5i + (52 - V{r) - h^v] 



a^JC^ - V2 - 2ia{JC,d, + JCydy) 



e — (5i + (52 



1 



e + 5i + S2 



(All) 



Applying the approximation given by Eq. (|A11[) to Eq. (jA10[) . we get from Eq. (|A10p the eigenvalue 
equation in the form: 



i-Si 4- 62 + V(r) + ^,,^y<^'^'~^'.-^M^f^+ JCydy) ^ ^^^^^2 P^JC^ - ^1 + 2^^^ + ICydy)\ ^ 
\ e-(5i-d2 e + (5i+(52 / 

= e^aa . (A12) 

Choosing the values for the coefficients a and /? to separate the coordinates of the center-of-mass (the 
wave vector IC) and the coordinates relative motion r in the Hamiltonian in the l.h.s. of Eq. (|A12p . we 
have 

e - (5i - ^2 
2^' 

/3 = '-^^ ■ (A13) 
Substitution of Eq. (Kn\i into Eq. (|JT2|) results in Eq. P^ . 
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Appendix B: Energy spectrum of an exciton 

Squaring both sides of Eq. , we get 

Fl = AN^-iFi . (Bl) 

Let us introduce the following notations: 

^1 = S1+S2 , 
v = 5i - S2 + Vo , 

Ci = 2-/N^{hvFf , (B2) 
which allows us to rewrite the Eq. (|Bip in the form: 



^--^^) =^^7^2- (B3) 
We can rewrite Eq. (|B3p as the form of the equation for e: 

+ Ae^ + Be^ + Ce^ + De + G = (B4) 

with the coefficients: 

A = -2v , 
B = v^ - {hvplCf - I? , 
C = (2//2 _ {hvFK.f )v - Ci , 

D = {{hvpJcY - v^)^? , 

G = {hvFK.fl-L^v . (B5) 

liv = {){5i = 52 & e = -Vq + e') then Eq. jBH has the form: 

(e2-(ftT;^/C)2)(e2-/i2)-Cie = (B6) 

with Ci < e(e^ — and /iufA^ ^ Z^- 

First, we assume Ci = and obtain eo = ±hvpK and eo = 

eo = ±hVFlC (B7) 
C = = (Ziy^^/C)^ + A , (B8) 

where A is very small correction. 

We substitute Eq. jBTj to Eq. (jB6P and neglect all the higher order quantities with respect to A: 



e = -Fo + ^{hvplCf + A2 ^ -T/o + \j -^^^^^ ^ (B9) 

which is not real, and, therefore, does not correspond to physical reality. 

Now let's consider the second solution of the zero-order expansion of Eq. ()B6|) : £0 = ±/i. Following 
the similar procedure, we obtain: 

(e^ - {hvFJCf){(^ - l-i^) - Cie = . (BIO) 
Now we substitute C = = A*^ + A that into Ea. (|B10[ ). where A is the small correction: 



-{hvFK.f)-Ci^/Ji? + A = Q , (Bll) 
and solve for A neglecting higher orders with respect to A: 

A = - - ^^)-^ - ^(M^ + i^VFKf) , (B12) 



which results in the fohowing energy dispersion: 



We expand this in powers of K. and obtain in second order 



/i2 
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Ci\ c 



\{hvFlCY) ■ (B13) 



e = -^o + A /^' + — ■ 1+ / X ■ (B14) 



This leads to Eq. 
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